Induced current in the presence of magnetic flux tube of small radius 
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The induced current density, corresponding to the massless Dirac equation in (2+1) dimensions 
in a magnetic flux tube of small radius is considered. This problem is important for graphene. In 
the case, when an electron can not penetrate the region of nonzero magnetic field, this current is the 
odd periodical function of the magnetic flux. If the region inside the magnetic tube is not forbidden 
for penetration of electron, the induced current is not a periodical function of the magnetic flux. 
However in the limit R — > 0, where R is the radius of magnetic flux tube, this function has the 
universal form which is independent of the magnetic field distribution inside the magnetic tube at 
fixed value of the magnetic flux. 
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I. INTRODUCTION 

The Aharonov-Bohm effect predicted in Refi has been 
investigated in numerous papers, see review^. Due to 
this quantum effect, the vector potential affects the elec- 
tron even in the case when the electron can not pen- 
etrate the region of non-zero magnetic field. This ef- 
fect has been investigated within non-relativistioi^— and 
rclativistio^ii^r— wave equations. In quantum field the- 
ory, similar effect has been investigated in Refs£~— . One 
of the consequences of the Aharonov-Bohm effect is scat- 
tering of a charged particle in the field of an infinites- 
imally thin solenoid. It was shown in Refsi^ that the 
Dirac Hamiltonian in the background of the field of in- 
finitesimally thin solenoid requires a self-adjoint exten- 
sion. This extension can be parametrized by one exten- 
sion parameter 9. As a result, the eigen functions of the 
Hamiltonian and the corresponding scattering cross sec- 
tion depend on this parameter. 

After recent fabrication of a monolayer graphite 
(graphene), RefJ^, new possibilities to study the 
Aharonov-Bohm effect have appeared. The single elec- 
tron dynamics in graphene is described by the massless 
Dirac equation in (2+1) dimensions, see RefSfiSr— . Us- 
ing the Green's function of this equation in the field of 
infinitesimally thin solenoid, the induced current density 
has been investigated analytically in RefJ^. Calculations 
have been performed exactly in a magnetic flux tube in 
the case when an electron can not penetrate the region of 
nonzero magnetic field. In the recent paper—, the local 
density of states in graphene was obtained in the field of 
infinitesimally thin solenoid. 

In the present paper, we derive the induced current 
density, corresponding to the massless Dirac equation in 
(2+1) dimensions, in the background of the magnetic 
flux tube of finite radius R. We investigate in detail 
the asymptotics R — > of this current density at fixed 
magnetic flux through the tube but for arbitrary distri- 
bution of a magnetic field inside the tube. We consider 
two cases. In the first case an electron can not penetrate 
the region of nonzero magnetic field, and in the second 
case this is allowed. Naturally, the result obtained in the 



first case is independent of the distribution of magnetic 
field inside the flux tube and coincides with the result of 
Refi^. More interestingly, in the second case the induced 
current density is also independent of the magnetic field 
distribution at R — > . In another words, we find the 
explicit value of the extension parameter 9, Refs^ i 11 ' 12 , 
and show that it is independent of the magnetic field 
distribution. 

The article is organized as follows. In Sec. |TT]we derive 
the Green's function of an electron in a magnetic field of 
the flux tube and obtain the induced current density at 
distances r ^> R. Calculation of the induced current 
density in the first and the second cases of the problem 
is performed in Sec. IIIII and Sec. IIV1 respectively. We 
conclude with a brief summary of our results. 



II. GENERAL DISCUSSION 

The induced current density, corresponding to the 
massless Dirac equation in (2+1) dimensions in a mag- 
netic field, has the form, RefJ^: 



3ind{r) = -ie 



{ |iTr{<TG(r,r| e )}. 
JC 2tt 



(1) 



where the Green's function G(r, r'\e) obeys the equation 

[e - cr ■ (p - eA(r))} G(r, r'\e) = 5{r - r') . (2) 

Here cr = (a x ,cr y ), and a% are the Pauli matrices; p — 
(j>x,Py) is the momentum operator, r = (x,y), A(r) is 
the vector potential of the magnetic field, e is the electron 
charge, the system of units h~c~ 1. According to the 
Feynman rules, the contour C of integration over e goes 
below the real axis in the left half plane and above the real 
axis in the right half plane of the complex e plane. The 
magnetic field of the flux is B(r) = B{r)v, where v is 
the unit vector directed along z axis, B(r) is a continuous 
function which obeys the condition B(r) = for r > R. 
It is convenient to choose the vector potential in the form 



A(r) 



<J>[v x r] 



V(r). 



(3) 
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where $ is the total magnetic flux through the tube, and Substituting the function D(r, r'\ie) in the form 
a function V(r) is 



V(r) 



2tt 



dyyB(y) . 



(4) 



Using the analytical properties of the Green's function, 
we deform the contour of integration over e in Eq. ([1} so 
that it coincides with the imaginary axis. Thus it is nec- 
essary to find the function G(r, r'\ie) which we represent 
as 

G(r, r'\ie) = [ie + a (p- eA(r))} D{r, r'\ie) , (5) 

where D(r,r'\ie) is the Green's function of the squared 
Dirac equation 



D(r,r'\ie) = 



and the ^-function as 



S(r-r r ) 



S(r - r') 



27tv rr 
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Jm(4>-4>') 



(7) 



(8) 



e 2 + (p - eA(r)f - ea z B(r) D(r, r'\ie) = -S(r - r') 



to Eq. (|6]), we obtain the following equations for the 
(6) functions T>m (r, r'\ie) 



dr 2 



ld_ 

r dr 



{™-iV{rjj l + 



jd(rV{r)) 
dr 



5{r - r') 



/ r r' 



(9) 



r 



Here s± = ±1, 7 = $/3>o, and $0 = 27r/e is the elemen- vanishes. As a result, the general solution of Eq. ([9]) in 
tary magnetic flux. In the region r, r' > R, the function this region has the form 
V(r) is V(r) = 1/r so that the term with s± in Eq. (|9]) 



Vg\rS\ie) = - (/| A |(|e|r<)i^| A |(|e|r>) + ^K w (\e\r)K w (\e\r')) , r,r' > R. 



(10) 
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Here I a {x) and K a (x) are the modified Bessel functions 
of the first and third kind, respectively, r> = max(r, r'), 
r< = min(r, r'), and A = m — 7. The coefficient in front 
of the first term in Eq. fj 10|) follows from the match- 
ing condition for the functions T>\n\r,r'\ie) at r = r' . 
The coefficients /3m depend on the distribution of the 



magnetic field in the tube. They can be found from the 
matching condition at r — R for the solutions (fT0|) and 
the solutions of Eq. © in the region r < R. 

Substituting Eqs. ^ and to Eq. (TT]), we arrive 
at the following representation for the induced current 
density at r > R 



3ind 



(r) = 



e]y x r] 

47T 2 7 



rn — — 00 



de 



~(vtiHr,r>\ie)+V(-Hry\ie)) - 1 (v(+)( r y\ie) -V(-\r,r'\ie) 



.(11) 
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Then we use the solution (flOl) for T>m \r,r'\ie) and finally obtain 
jind{r) = jln d (r) +3-n d (r), 

JiniM = "^2^ E / de/| A |(|e|r)if| A |(|e|r) = ( - - | 7 | tan(7r 7 ) . 



m— — oo 



jind( r ) 



e[i/ x rj 
27r 2 r 



E 



rn— — oo 



rfr 



(/3(+)X| A | +AI (| e |r) -/?(-)^| A |_ A1 (|£|r)) %|(|e|r) : 



r 



(12) 
(13) 

(14) 



where /i = sign(A), the quantity 7 is the fractional part 
of 7, | 7 | < 1. It is defined as 7 = 7 — n for 7 > 0, and 
7 — 7 + n for 7 < 0, n is the maximal integer number less 
than | 7 |. The first contribution, j^ d (r), coincides with 
the result of RefJ^ for the induced current density in 
the field of infinitesimally thin solenoid. This result has 
been obtained under assumption that an electron can not 
penetrate the region of nonzero magnetic field. The con- 
tribution is a periodical function of the magnetic 
flux, it depends on the fractional part of the quantity 
7. It is equal to zero at the integer values of 7 but also 

at half-integer values of this quantity. The contribution 

(2") 

3ind( r ) depends on the radius R and distribution of the 
magnetic field inside the magnetic tube through the co- 
efficients Pm^- These coefficients are the functions of e, 
R, 7, and m. 



The corresponding boundary condition is ji n d(R) = 0. 
In this case T>m\r, r'\ie) = T>m\r, r'\ie) and we obtain 
from Eq. (II 1 p for the induced current density 



jind{r) 



X deVM{r,r'\ie). (15) 
, Jo 



The boundary condition gives the following form of the 
coefficient /3^ } in Eq. (HUJ) 

= -I\x\{x)/K w {x), x = \e\R. (16) 



At X -C 1, the coefficient j3m in the leading order in x 
reads 



III. THE CASE OF PURE AHARONOV-BOHM 
EFFECT 



w ~ 21 " 2|A| I A I #\> 

~~ I*(|A| + 1) 



(17) 



In the pure Aharonov-Bohm effect, an electron can not 
penetrate the region r < R of nonzero magnetic field. 



Using this expression, we obtain the asymptotics of the 

? -(2) / 
find v 



term j ind {r) in Eq. (fT4]) at large distances r » R 



? -(2) 



(r) 



e[^xr] ~ A|A| £(|A| + l/2)r(2|A| + l/2) ^i?N 2X 



27rr 3 



rn— — 00 



r 3 (|A| + i) 
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The contribution j ind {r) is equal to zero at integer and half-integer values of 7 due to cancelation between different 
terms. For other values of 7 we introduce the quantity A — rh — 7, where fh minimizes the value | 7 — m|, and obtain 
jind( r ) m the leading and next-to-leading order in R/r 



3ind{r) 



e[v X r] 



1 - m)' ta„(»fl + 2A|X| r(|A| + l/2)r(2|A| + l/2) / A\ 2W 

2 "7 v ,; 11 r 3 (|A| + l) \2r 



(19) 



The leading term j (1) d (r) tends to zero at r -> 00 as 1/r 2 , IV. THE CASE WHEN AN ELECTRON CAN 
and the correction (r) diminishes as (5/r)^l/r». PENETRAT mIgNE R T E IC I ?iel°d N ° NZER ° 



In this section we assume that an electron can pene- 
trate into the magnetic tube. Let us start with a sim- 
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pie case of the magnetic field distribution of the form 
B(r) — B 8(R — r), where 6{x) is the step-function. In 
this case the function V(r) reads 

V(r) = ^-6{R-r) + - 6(r - R) . (20) 
R z r 

This function is continuous one at r = R so that the 
matching conditions for the functions T>m (r, r'\ie) at r = 
R, following from Eq. are 

V£\R - 0, r'\ie) = V£\R + 0,r'\ie) , 



&D£\ry\ie) 



Or 



dV^\ry\ie) 



r=R-0 



dr 



(21) 



r=R+0 



The solution of Eq. (0) in the region r < R, r' > R has 
the form: 



V£\r,r'\ie) 



C^VK^'exp 



,(±) 



x 1 F 1 

e 2 R 2 
4| 7 | 



,(±) 



, \m\ + 1, 



+ 



\l\r 2 
2R 2 

\l\r 2 
R 2 

■ym + 7 1 



(22) 



s±7 
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where C m (r') are some functions which can be found 
from the conditions (j2"Tj) . iF% (a, b, x) is the confluent Hy- 
pergeometric function. Using Eqs. (fTOl) . (|2T1) . and (|22|) 
we obtain for the coefficients /3™ 

(±) = [6J| A |(aO + xl w+1 (x)] - ^ri^jx) 

F[bK w {x)-xK w+l {x)]-2\~f\FK w {x) ' 

^ = 1J F 1 (a( ±) ,|m| + l,| 7 | N 



m 



&=|A| + | 7 |-|m| 



i,H 

x = leli?. 



2,|7l 



(23) 



It follows from Eqs. (1231) and (TT4")) that the induced cur- 
rent density is the odd function of the parameter 7, which 
is a consequence of the Furry theorem. Below for sim- 
plicity we assume that 7 > 0. 

Using Eqs. (f2"3")l and (TT4")) it is easy to obtain the asymp- 
(2) 

totics of the term j\ nd {r) at distances r >• R. In this 
case, the main contribution to the integral over e in Eq. 
(TT4")) is given by the region e ~ 1/r so that a; = |e|i2 -C 1. 
The main contribution to the sum over m in Eq. (TTJ 

?(+) 



is given by the term proportional to (3m 1 with m = [7], 

this 

at x <C 1 is 
■ sin(7T7) . 



where [y] means the integer part of y. For this m the 
asymtotics of the coefficient /3„ 

2 

7T 



(+) 



(24) 



Substituting this expression to Eq. (fT4"f and taking the 
integral over e we find 



Airr 3 



|7| tan(7T7) . (25) 




FIG. 1: Dependence of the function F( 7 ), Eq. (|26J) . on the 
7 = e$/(27rc), where $ is the magnetic flux. 



Taking a sum of j\n d {r) and J,-^( r ) we nnau y obtain 



jind(r) 



e[v x r] 



F( 7 ). 



167rr 3 

F( 7 ) = (1 - 47 2 ) tan(7T7) 



(26) 



The function ^(7) is shown in FigjTJ 

As it was shown in Ref J^, the induced current in the 
case of pure Aharonov-Bohm effect is the periodical func- 
tion of a magnetic flux. The reason is very simple. When 
an electron can not penetrate the region of nonzero mag- 
netic field, the integer part of 7 can be eliminated from 
the wave equation by the gauge transformation. In oppo- 
site case considered in this section, we can not eliminate 
the integer part of 7. As a result, as it is seen from 
Eq. (|25|) and FigfTJ the induced current density is not a 
periodical function of the magnetic flux. 

The result (|26|) . obtained for the magnetic field B(r) = 
Bq8(R, — r), coincides with the corresponding result of 
Refiii for the extension parameter — ir/2. Now we 
show that the leading asymptotics of the induced current 
density at R/r -C 1 is independent of the distribution 
of magnetic field at the fixed magnetic flux. As it was 
pointed out above, at R/r 1 the main contribution 
to the integral over e in Eq. fT4")) is given by the region 
e ~ 1/r so that |e|i? <C 1. For r,r' > R the solution 
(fTU)) is valid for any distribution of magnetic field. We 
use this solution for r' > r > R. For r' > R > r we can 
neglect the term e 2 in Eq. © and write the equations 

for the functions T>m(r, r'\ie) in the form: 

A m B m V^(r,r'\ie) = 0, 

B m _ 1 A m _ 1 v(-\r,r'\ie)=0, 
d m + 1 



Am — tt~ 

or r 



7 V(r), 



dm 

B m =x + lV{r) 

or r 



(27) 



The regular at r = solutions of the Eq. (|2"T)l have the 
form 
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W (r, r'lie) = #> (r>« exp {- 7 £ Vfo)*} { }^ ^ (2ro+1) exp ^ /; & ^ ° ; (28) 

i?k- ) (ryiie) = cir ) (r')r-«p{ 7 /V(^) { f * 9m_1 ex p {- 2 ^^ n#} , m>o, (29) 



Here £m (r') are some functions of r'. Substituting the 
solutions (fTU)) . (|2"51) . and ([25)) to the matching condition 
(f2"T|) and using the asymptotics of the functions I\(x) 

and K\(x) at small a; we obtain the coefficients j3^\ In 
the limit \e\R — > 0, only /3,*t' does not vanish and co- 
incides with the result obtained for the magnetic field 
B(r) = B 9(R - r), see Eq. 423}. Thus the leading 
asymptotics of the induced current density is indepen- 
dent of the distribution of the magnetic field in the mag- 
netic flux tube. In our consideration we do not use any 
assumptions concerning the value of magnetic flux. 

V. CONCLUSION 

We have investigated the induced current density, cor- 
responding to the massless Dirac equation in (2+1) di- 
mensions in a magnetic flux tube of small radius. This 
problem is important for graphene for which all formu- 
las for the induced current density obtained in our pa- 



per should be multiplied by the factor Avf, the Fermi 
velocity is Vp sa 10 6 m/s rj c/300. In the case of pure 
Aharonov-Bohm effect, when an electron can not pene- 
trate the region of nonzero magnetic field, this current is 
the odd periodical function of the magnetic flux. In the 
opposite case, when the region inside the magnetic tube 
is not forbidden for penetration of electron, the induced 
current is not a periodical function of the magnetic flux. 
However in the limit R — ¥ this function has the uni- 
versal form which is independent of the magnetic field 
distribution inside the magnetic tube at fixed value of 
the magnetic flux. Thus it is shown that the extension 
parameter 9 introduced in the previous papers has the 
universal value tt/2. In the same way it is easy to show 
that the wave functions corresponding to the massless 
Dirac equation in (2+1) dimensions in a magnetic flux 
tube have also the universal form in the limit R — > 0. 
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